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ABSTRACT 

With the help of the Laplace-Lagrange solution of the secular perturbation 
theory in a double-planet system, we study the occurrence and the stability of 
apsidal secular resonance between the two planets. The explicit criteria to pre- 
dict whether two planets are in apsidal resonance is derived, which shows the 
occurrence of the apsidal resonance depends only on the mass ratio (mi/m2), 
semi-major axis ratio(ai/a2), initial eccentricity ratio (eio/e2o), and the initial 
relative apsidal longitude {w2q — Wio) between the two planets. The probability 
of two planets falling in apsidal resonance is given in the initial element space. 
We verify the criteria with numerical integrations for the HD12661 system, and 
find they give good predictions except at the boundary of the criteria or when 
the planet eccentricities are too large. The nonlinear stability of the two planets 
in HD 12661 system are studied by calculating the Lyapunov exponents of their 
orbits in a general three-body model. We find that, two planets in large eccentric- 
ity orbits could be stable only when they are in aligned apsidal resonance. When 
the planets are migrated under the planet-disk interactions, for more than half 
of the studied cases, the configurations of the apsidal resonances are preserved. 
We find the two planets of HD12661 system could be in aligned resonance thus 
more stable provided they have — ~ 180°. The applications of the criteria 
to the other multiple planetary systems are discussed. 



Subject headings: celestial mechanics-method: analytical and numerical -planetary 
systems-stars: individual (HD12661, 47 Ursae Majoris) 
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1. INTRODUCTION 

The detecting of extrasolar planetary systems reveals fruitful results during the past 
years. More than 100 extrasolar planets has been inferred by the Doppler radial velocity 
measurements to the solar-type stars (California and Carnegie Planet Search 2003), among 
them 10 multiple-planet systems are confirmed. For a multiple-planet system, the dynamical 
stability of the system under planetary interaction is an important issue concerning the 
dynamical evolution as well as the possible existence of habitable zone of the system. 

There are many effects which can affect the stability of a multiple-planet system. For 
the orbits of planets with small or modest eccentricities and inclinations, mean motion 
resonances between planets can sometimes lead to stable configurations. Another effect is 
the secular resonance between the planets. An apsidal resonance occurs when the relative 
apsidal longitudes of the two orbits Atu librates about (aligned resonance) and n (anti- 
aligned resonance) during the evolution. Due to the aligned apsidal resonance, the two 
planets on elliptic orbits can greatly reduce the possibility of close encounters, thus it is 
believed that aligned apsidal resonance can stabilize the interacting planets. 

For the 10 multiple planetary systems observed to date (GJ876, 47UMa, HD82943, 
HD12661, HD168443, HD37124, HD38529, HD74156, Upsilon Andromedae, and 55 Cancri), 
the best fit orbital parameters inferred from the radial velocity observations show 6 pairs 
of planets could be in apsidal resonance: HD82943(Gozdziewski & Maciejewski 2001), Up- 
silon Andromedae c and d(Chiang, Tabachnik, & Tremaine 2001), GJ876(Lee & Peale 2002), 
47UMa(Laughhn, Chambers, & Fischer 2002), HD12661(Gozdziewski & Maciejewski 2003; 
Lee & Peale 2003), and 55 Cancri b and c(Ji et al. 2003). Ubiquitous as it is, the apsidal 
resonance phenomenon is worthy to be studied in detail. In this paper, we are interested 
when the apsidal resonance occurs and whether it really leads to a stable configuration be- 
tween planets, since an anti-aligned resonance could lead to close encounters between orbits 
with large eccentricities . 

For the occurrence of apsidal resonance, Laughlin et al. (2002) gives a criterion based 
on the Laplace-Lagrange solution of the secular perturbation system. In this paper, the 
criterion is represented to a more explicit form in section II. The probability that the apsidal 
resonance happens is also derived according to the criteria. In section III, we study the 
stability of orbits in apsidal resonances by calculating the largest Lyapunov exponents of 
the orbits with a general three-body model. The behavior of orbits in migration are studied 
with a torqued three-body model. The conclusions and the applications of the criteria to 
the other multiple planetary systems are discussed in the final section. 
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2. LOCATIONS OF APSIDAL RESONANCE 

In this section we derive the exphcit criteria under which the apsidal secular resonance 
may occur. To that aim, the hnear secular perturbation theory is employed for two in- 
teracting planets under the attraction of the host star. For a planetary system with two 
planets, hereafter we denote all the quantities of the host star, the inner and outer planets 
with subscripts "0", "1" and "2", respectively. So the three bodies have masses mo, mi, m2, 
respectively, where mi,m2 <^ mo- In the present study we address the coplanar problem 
only, so the two planets are on the orbits with osculating orbital elements {ai, ei. zui, Mi) 
and {a2, 62, M2), respectively, where a, e, w, M are the semi-major axis, eccentricity, lon- 
gitude of pericenter and mean anomaly of the orbit, respectively. We adopt the commonly 
used unit system, i.e., the mass unit is the solar mass, the length unit is lAU, and the time 
unit is lyr/(27r). 
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2.1. Linear seculeir perturbation theory revisited 

We start with the linear secular theory following Murray and Dermott (1999). For the 
coplanar case, the disturbing function for the motions of planets mi and m2 are given as: 

Ri = riial l^Aiiel + ^126162 cos(roi - rxj2)] 
R2 = n2al [\A22el + ^216162 cos(g7i - ^2)] 

where ni,n2 are the mean motion of planets mi and m2, respectively; Aij are elements of 
matrix given by 

^ Cl -coci \ ^2) 

-C0C2 C2 / 

where Cfe > (/c = 0, 1, 2) are functions of oi, 02, mo, mi, m2 defined as 

co= hf{a)/h^l\a)^\a{l-\a') 

ci= i^^T^^^i^ (3) 
\n2^ahf{a) 

with h^^ {a){i = 1,2) being the Laplace coefficients, and a = 01/02 < 1. Fig.l shows the 
2 

approximation of Cq by the above formula. Quantitatively, the error of approximation is less 
than 1% for a < 0.66, and less than 5% for a < 0.91. So the approximation is quite well for 
the study of the planetary system. Moreover, we define 

C2^ ln2mi(mo + mi)^^^,,2^ (4) 
Cl a nim,2{mo + 1712) 
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with q — mi/m2 and the terms with orders of 0{^),0{'^) or higher are neglected in the 

mi m2 
mo ' mo 



above approximation, since ^ ~ 



10 ^ in the planetary systems. Denote 5'i,g'2 as the 

cos 6i 



two eigenvalues of matrix (2), and the corresponding eigenvectors are S,, 



Oi e (-7r/2,7r/2) and 
cos 9i 



CqCi 
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Define 



S-i = |[(ci + ca) + a/(ci - C2)2 + 4c^ciC2' 



sin 



1,2) 



, where 



(5) 



5'2 = |[(Cl + C2) - Vl^l - C2)^ +4cgciC2]. 

= 2^(l-e+A/(l-0^ + 4cgO- 



P2 



tan 6*9 



(6) 



(7) 



So pi < and p2 > 0, or — | < 6'i < < ^2 < f • The scaling factor Si (i=l,2) can be 
expressed in terms of initial eccentricities eio, 620 and Awq = W20 — zuiq: 



[/^2^io~^P2eioe2o COS Aroo+e|o]^/^ _ F 

|pi-p2|cos9i |pi-p2| cosSi ' 

[plelo-^Pieioe20 cos Aroo+e|o]^/^ _ G 

\pl-p2\c0sB2 |pi -P2 1 008 6*2 ' 



(8) 



The secular system with disturbing functions (1) is integrable and the solutions can be 
written as: 



ei 

62 

6162 sin Aw 
6162 cos Azu — 



^[F2 + 2FGcosA^ + G2]V2 
^[p?F2 + 2pip2FGcos AV^ + p2G2]V^ 



(9) 



P1-P2 
1 



[piF2 + (pi + p2)FG cos A^ + P2G2]_ 



(Pl-P2)2| 

where A'^ — ip2 — i^i — {g2t -\- P2) — {git + with t the time and /32 given by: 

sin/3i = j^{hiop2 - /i2o), sin/^s = -^(/iioPi - ^20), 

cos/?i = j^{kioP2 - k2o), cos/52 = -^{hopi - /C2o)- 



(10) 



where hio — eiosinzuio, kio — Ciocoszuio, {i — 1,2). From (9), it's easy to verify that the 
evolution of Ci, 62 obeys an integral: 
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where D is a constant which depends only on the initial parameters. Moreover, from (8) and 
(9), the maximum of ei and minimum of 62 occur at cos A'^ = 1 (as pi < 0), with values: 

. _Z + G_ ^ _ \piF + p2G\ 

The minimum of Ci and maximum of 62 occur at cos = — 1, with 

_ \F-G\ _ p^G - p,F 

^Imm I I 5 ^2max i i ■ v / 

\Pl-P2\ IPI-P2I 

Thus, we can obtain the maximum excursions of ei and 62 for any given eio, 620, Aotq as 
follows: 

+ _ {P2G - p,F) - \p,F + p,G\ ^^^^ 

IP1-P2I ' IP1-P2I 



2.2. The explicit criteria 

With the help of the last equation of (9), the criterion in Laughhn et al. (2002) for the 
apsidal resonance can be expressed as 



{p^+P2)FG 



p,F^ + P2G^ 



< 1 (15) 



Since when S < 1, the values of Azu can not reach 7r/2 or 37r/2 (thus cosAtu ^ Q), so it 
must librate about or tt. On the contrary, when 5* > 1, it is possible that Azu will reach 
71 / 2 or 37r/2 , thus it will circulate in [0, 27r]. 

The equation (15) equivalent to, after some algebra manipulations: 

F F 

— > max(l, -P1/P2), or, < — < min(l, -P1/P2), (16) 
In view of (8), the above relations are equivalent to : 

— < cosAoto, (17) 

eio Pi + P2 

or, 

'J^> P^+P^ \ >o. (18) 
eio 2 cos Azuq 

By substituting (3) (4) and (7) into the above expressions, we finally obtain: 

620^ 5 ga^/^(l-X) 

— < ~o — 1 — cosAtiJo, 19) 

eio 2 1 — go;^/^ 
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or, 



620 2 1 



> 0. 



(20) 



Cio 5 Q!(l — |q;2) cos Aot( 



'0 



These are the cxphcit criteria for the occurrence of the apsidal secular resonance. Equations 
(17) and (18) are obtained with the hnear secular perturbation theory, while to get (19) and 
(20), we use the approximations of cq and ^ in (3) and (4). 

We call the libration region defined in (19) the down-libration region, and that defined 
in (20) the up-libration region. Whether the down-libration or up-libration is the aligned 
or anti-aligned resonance depends on the sign of (1 — qa^^"^). For qa^^"^ < 1, down-libration 
occurs only when 7i/2 < Azuq < 37r/2 (thus it is the anti-aligned resonance) and up-libration 
occurs when Awq < 7r/2, Awq > 37r/2 (the aligned resonance). This is the case of HD12661 
system. If ga^/^ > 1, the conclusions are reversed, which is the case of 47UMa system. For 
the critical case qa^^"^ — 1, according to (19) and (20), all the orbits are in libration except 
those with Aujq — 7r/2,37r/2. Fig. 2 shows a typical phase space and the evolution of an 
orbit in the linear secular perturbation system (1). The jumps of Azu are due to the cross of 
origin in the (cjCostUj, ejsintUi) plane. Fig.3 shows the resonance region in the 620 — Azuq 
plane defined by (19) (20) with different Ciq. The parameters a and q are taken from the 
planetary systems HD12661 and 47 UMa(listed in Table 1-2). The boxes around the present 
configuration dots show the uncertainties of the elements (also listed in Table 1-2). 

There are two limit cases for the criteria (19) and (20): 

(1) q; — > 0. The minimum 620 for up-libration and maximum 620 for down-libration can 
be obtained by setting cosAciJ = 1 in (19) and (20). When a — > 0, the minimum 620 for 
up-libration tends to very large and maximum 620 for down-libration tends to zero. So when 
the two planets are far away, the libration regions in 620 — Awo plane can be negligible small. 

(2) Cio — > 0. This happens when one of the planets is in a near-circular orbit, and 
it is just the case discussed in Malhotra (2002), namely for two planets initially in nearly 
circular orbits, an impulse perturbation may impart a finite eccentricity to one planet's 
orbit. When Cio — > 0, criterion (20) is always fulfilled if —n/2 < Awq < 7r/2 for qa^^^ < 1 
or 7r/2 < Azuq < Sir/ 2 for qa^^"^ > 1. The boundary curves are the limits of those with 
cio tends to zero in Fig. 3a or Fig. 3b, with both the minimum 620 for up-libration and the 
maximum 620 for down-libration tend to zero. Thus half of the e2o — Azuo plane is the possible 
resonance region with up-libration. So the probability of these two planets captured into 
apsidal resonance tends to 50%. This conclusion agrees with that of Malhotra (2002). It's 
possible that these two planets are captured in anti-aligned apsidal resonance, which depends 
on the sign of 1 — qa^^'^. Similarly, when 620 — 0, criterion (19) is always fulfilled, thus half 
of 620 — Awo plane are possible resonance region with down-libration, and the probability of 
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the two planets in apsidal resonance (either aligned or anti-aligned resonance) tends to 50%. 

To compare the above libration regions obtained by the criteria with those calculated 
from the secular perturbation system, we integrate the orbits of the secular perturbation 
system for the HD12661 system. Fig.4 shows the diagrams of orbits in 62 — Aro plane, the 
initial values of the studied orbits have eio = 0.1, with /S^voq — for Fig.4a and /\.vdq — tt 
for Fig.4b. According to the criteria(19) and (20), 620 > 0.022 for the aligned resonance at 
Agjq = 0, and 620 < 0.375 for the anti-aligned resonance at A.vjq — tt, which coincide with 
those given from the secular perturbation system. 

Fig. 5 shows the contours of the excursions of 61,62 for different 620 and Aotq. As we can 
see, generally the orbits in apsidal resonance have relatively smaller A61, A62, especially the 
turning points of the contour curves lie in one of libration boundary curves. Thus from the 
linear secular perturbation theory, the orbits in apsidal resonance, either in aligned resonance 
or anti-aligned resonance, are more stable than those in non-resonance region. 



2.3. Area of the libration region 

From the above criteria, we can calculate the probability that the two planets fall in 
apsidal resonance in the space of initial orbital elements . We define the probability as the 
area of the libration region in the 620 — At^o plane for a given Ciq. For the down-libration 
case, according to (19), it's possible that the peak value of the 620 — Aotq curve can be above 
unit for larger 610 (as the 610 = 0.35 case in Fig.3a, and the 610 = 0.50 case in Fig. 3b). We 
set /^VDcL as the half width of the down-libration region where the boundary curve reaches 
620 = 1- Fig. 3a and Fig.3b show Aro^ for the 610 = 0.35 and 610 = 0.50 curves, respectively. 
Define 

_5 eioga^V^(l-X) 
^'-2 |l-gaV2| ' (21) 

then, 

A = / arccos(l/Qd), if > 1 /22l 
I 0, if Qd < 1 ^ ' 

and the area ratio of the down-libration area to the total area of 620 — Atoq plane is. 

Pa = li^wa + Xfr:^^"^"^ Q.cos A^ dl^w] ^23) 
= ^[^Wd + (9d(l - sin AtUd)] 

where the lower integration limit is the beginning point of the down-libration region in Aotq- 
axis {/^Wa = f for qa^/'^ < 1 and Aro^ = ^ for qa^/'^ > 1). Fig.3a and Fig.3b show the 
/^VDa for the 610 = 0.35 and 610 = 0.50 curves, respectively. As one can see, the area ratio 
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for the down-libration increases linearly with eio when eio is small, since in the interval one 
has = in equation(23). However, for larger cio, the increase of area ratio is no longer 
linear since Aro^ ^ 0. 

Similarly, for the up-libration resonance, if we set Aci7„ as the half wide of the up- 
libration region when the boundary curve meets 620 = l(see FigS) and define 

5 a[l — ^a^) 

so, 

Ati7„ = arccos((5u), (25) 
and the area ratio of the up-libration area to the total area of 620 — Atuo plane is. 



= 1[A^„-Q„ln(ii|^)] 



(26) 



where the lower integration limit Aro^ is the center of the up-libration region (Aro^ = for 
ga^/^ < 1 and Aro^ = tt for qa^^'^ > 1, see Fig.3). 

In the early evolution of planetary systems, both q and a may vary due to planetary 
formation and migration. We fix eio — 0-35 for the HD12661 systems, and see the variation 
of libration area ratio with a or q. Fig. 6 shows the variation of libration area ratios with 
a and q. In Fig. 6a q^ 1.46 is fixed as the observed value. One can see the ratios increase 
with a before they reach the maximum (unit) at a„iax = Q""^ ~ 0.47, which is the critical 
case, and then decrease. In Fig.6b a ~ 0.32 is fixed as the observed values. Again the curves 
reach the maximum at g = o;"^/^ ~ 1.76, the critical case. 



3. STABILITY OR ORBITS IN RESONANCE 

Since the linear secular perturbation theory is an approximation to the real three-body 
system, the above criteria obtained from the linear perturbation theory has its limitation. 
To apply the linear criteria to the predicting of the apsidal secular resonance, we integrate 
the orbits in a general three-body (co-planar) system, where the longitudes of the ascending 
nodes and the inclinations of the two planet orbits are assumed to be zero {Qi = ^2 = 0, 
h = ^2 = 0) in the paper. We adopt the RKF7(8) (Runge-Kutta-Fehlberg) integrator with 
adaptive step-sizes to integrate the orbits. Generally the step is alternated between 0.00625yr 
and 0.0125yr, so there are 80-160 steps in a period of planet orbit with a semi-major axis 
of lAU, and the final error of the Hamiltonian of the three-body system after lOM years' 
evolution is less than 10~^. 
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We define an index to indicate whether an orbit is in hbration region or not. Choose 
a serial of discrete time during the evolution of orbits (for example, every 12.5 years), give 
an index J„ for each time so that J„ = if at that time — 7r/2 < AciJ < 7r/2 and /„ = 1 if 
7r/2 < A-a7 < 37r/2. Then the average values of J„ over very large n, denoted by < /„ >, 
shows roughly the character of the orbit during the studied period of time according to 



Fig. 7 shows the index for the HD12661 planet system for 10 x 51 orbits in the interval 
[0,0.5] X [0, tt] of the 620 — Azuq planc(the same for the following calculations), the initial 
eccentricity is Cio = 0.10 in Fig.7a and Cio = 0.35 in Fig. 7b. As one can see, most of the 
orbits in the libration region predicted by linear criteria (19) and (20) arc in real libration 
region when 620 is small. The discrepancies between the linear system and the three-body one 
mainly occur on larger 620 and the boundary between the libration and circulation region. 

Next we want to see whether the orbits, cither in apsidal resonance or not, may have 
different stability in the general co-planar three-body model. The stability of an orbit in a 
Hamiltonian system is related with the topology (regular or chaotic) of the phase space, so 
we calculate the largest Lyapunov Characteristic Exponent (LCE) to indicate whether the 
corresponding orbit is in a regular or chaotic region. The LCE at finite time x(t) is calculated 
for few orbits up to i = lOMyr (denoted as xy), and for most orbits up to t = IMyr (denoted 
by Xe)- Fig-8 shows the LCEs of four orbits in the four different kinds of region in the phase 
space. For curve (a), x(^) decrease linearly with t, thus the orbit corresponding to curve 
(a) has zero LCE, and is in a regular region. Curve (b) shows a very small but non-zero 
LCE, so the orbit corresponding to (b) is in a very weak chaotic region. Both curves (a) 
and (b) have Xe ~ 10~^yr~^, x? ~ 10~^yr~^. Curve (c) tends to a constant value, with 
Xe ~ x? ~ lO^'^yr^^, so the orbit corresponding to (c) should be in a strong chaotic region. 
The orbit corresponding to curve (d) is unstable with the outer planet escape before lO^Myr, 
and x(^) for such an orbit is generally great than 10~^yr~^ before escape. Thus by calculating 
Xe (unit: yr~^), we can tell at least three different kinds of orbits: 




1 

others 



aligned libration 

circulation 
anti — aligned libration 
mixed 



(27) 




regular or weak chaos 
strong chaos 
unstable 



(28) 



We calculate xe for the HD12661 system with initial eccentricity eio = 0.10 in one run 
and Cio = 0.35 in another run, and the other initial parameters are taken as the observed 
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values. Fig.9 show the results. The boundary curves for the corresponding eio are also 
plotted in the diagram. We find for small eio, the orbits, whether in aligned resonance, 
anti-aligned resonance or non-resonance regions, do not show much difference about the 
LCEs. So in this case, whether an orbit is in apsidal resonance do not affect much on its 
stability. However, for larger initial eio, orbits in aligned resonance seem to be more stable 
since they have much lower LCE as compared with those in the anti-aligned resonance or 
circulation regions with same 620- This example shows for larger cio, 620, planets in aligned 
resonance regions would be relatively more stable. Prom Fig.9b, we can also see the present 
configuration of HD12661b and HD12661c is in the boundary of a chaotic region, if Aotq is 
set to 130.7°, which is symmetric with the observed value —130.7'^. This conclusion has been 
obtained by Kiseleva-Eggleton et al. (2002). 

Finally, we study the role of apsidal resonance on the stability of the orbits when the 
planets are migrating. In the early stage of planet evolution, the protoplanets and the 
stellar-disk might be coexisting and interacting, thus planet migration might happen due to 
nebular tides(see, e.g.. Ward 1997). Wc adopt the torqued three-body model as in Laughlin 
et al. (2002), and for the sake of simplicity, we consider the case that only the outer planet 
experiences an azimuthal torque due to the planet-disk interaction. We take the azimuthal 
acceleration as /2 = —2 x 10~^AU^yr~^ (which is smaller than that used in Laughlin et al. 
(2002), because here we concern the quahtative evolutions only), and study both the forward 
{t > 0) and backward {t < 0) evolutions of orbits under this acceleration. We calculate the 
orbits for the HD12661 system with initial Cio = 0.35, and Mio, M20 are random chosen. All 
the other initial parameters of orbits are taken from table 1. The evolution time span is 50000 
years. We find in all the studies cases(both forward and backward), the semi-major axis of 
the inner planet do not have secular changes, just as the results in Laughlin et al. (2002). 
Fig. 10a shows the index for the orbits during the evolution and Fig. 10b shows the 02 at the 
final time for the forward case. Some of the orbits which were initially in apsidal resonance 
region become mixed, though, according to Fig. 6a, the libration regions are enlarged due to 
the increase of a (from 0.32 to approximately 0.41) by planet migration. For orbits in the 
aligned resonance, they are more stable and with smaller final 02, while for orbits in anti- 
aligned resonance with larger 620, or in circulation region, they tend to be more unstable 
due to close encounters and thus have larger final 02, and most of them will escape soon 
in the following evolutions. Fig. 11 shows the backward case. Wc find the conclusions are 
more or less similar. Though the libration regions shrink due to the decrease of a by the 
planet migration (from 0.32 to approximately 0.23) in this case, the configurations of apsidal 
resonance are generally preserved during the migration. Thus most of planet systems which 
are observed in apsidal resonance now may be in resonance before the migration begins. 
Again, the orbits in aligned resonance seem to be more stable in the sense of having modest 
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final a2, which may be the results of less close encounters between the two planets. 

4. CONCLUSIONS AND DISCUSSIONS 

In this paper, we have studied the occurrence as well as the stability of the apsidal 
resonance. The apsidal resonance occurs when the equations (19) (20) are fulfilled. We find 
the occurrence of apsidal resonance depends only on the mass ratio q — mi/m2, semi- major 
axis ratio a — ai/a2(in secular systems, ai, a2 are constants), initial eccentricity ratio eio/e2o 
and the relative apsidal longitude W2o—'^w of the two planets. The criteria are based on the 
Laplace-Lagrange secular solution of linear perturbation theory. Based on these criteria, the 
ratio of librating to non-librating orbits in the 620 — Atug plane can be obtained analytically, 
which is given in (23) (26). We also find for two planets on the orbits with large eccentricities, 
they can be in a stable configuration only when they are in aligned apsidal resonance. When 
the planets are migrated under the planet-disk interactions, more than half of the studied 
orbits preserves the configurations of apsidal resonance. 

The linear secular perturbation theory is applicable only when the two planets are not 
in a lower order mean motion resonance. Since in lower order resonances, the variations of 
zui and ZU2 are not guided by the secular dynamics, but the resonance angles. For example, 
for the j : {k — j) resonance, if the two resonance angles 9i — jXi + {k — j)X — kzui, 
62 — jAi -\- {k — j)\ — kvD2 librate around or tt, then the relative longitude of pericenter 
Aro = (^1 — 62) /k must librate around either or tt. Thus we think in this case, the 
lower order mean motion resonance and the apsidal resonance are not independent, and the 
former one guide the dynamics. In the case that only 9i librates, we believe that the apsidal 
resonance is very difficult to occur, since in this case Wi is guided by 6*1, thus it can not have 
similar variations with W2- 

Beauge et al. (2002) find that, there may exist some asymmetric stationary solutions in 
mean motion resonance region, where both the resonant angles and Atu are constants with 
values different from or tt. We think such kind of solutions are due to the mean motion 
resonance and can only exist in the resonance regions, since such kinds of apsidal resonance 
solutions with Aro librates about constants with values different from or tt can not be 
found in the linear secular perturbation system. 

For the ten known multiple-planet systems (see Table 8 of Fischer et al. (2003) for 
a list of the elements), the situation of whether apsidal resonance happens between their 
planets can be classified roughly into three groups (see Table 3 for the extensions of Aotq 
when apsidal resonance would occur for the observed a, ^): 



- 12 - 



(1) planets both in apsidal resonance and mean motion resonance. 55 Cancri b and c 
(Fig. 12a), GJ876 b and c , HD82943 b and c are in apsidal resonance since they are in mean 
resonances 3:1, 2:1 and 2:1 respectively. In fact, HD82943 b and c can in apsidal resonances 
without in the mean motion resonance, while GJ876 b and c are near the boundary of 
libration according to the linear secular dynamics(Fig.l3a). 

(2) planets in apsidal resonances far away from lower order mean motion resonances. 
HD 12661 b and c, 47 UMa b and c, Ups And c and d(Fig.l2b) are in this type. They are in 
apsidal resonance without the existence of any strong mean resonances. Moreover, HD12661 
b and c seems to be in anti-aligned apsidal resonance, which is in the boundary of a chaotic 
region. 

(3) planets not in apsidal resonance either due to the negligible small libration region 
in the 620 — Awq plane or without suitable Awq. The two planets in HD38529, HD168443, 
HD74156 are not in apsidal resonance, since the libration regions in the 620 — Azuo are 
neghgible small for the observed q,a and 62/61 (Table 3). HD37124 b and c have large 
aligned libration area for the present parameters, and their eccentricities are not small, but 
they are not in resonance due to the present values of Azuq if Jl^ = = are assumed 
(Fig.l3b). 

However, due to the unknown of the inchnations and the longitude of ascending nodes 
in the orbital fit from the observation data, it is still too early to make conclusions for 
some planet systems whether the planets are in apsidal resonances. For example, for the 
HD12661 b and c, they are believed to be in the anti-aligned libration now. This is achieved 
by assuming Qi = Q2 = 0, and they are in the boundary of a chaotic region. Alternative 
choices of Qi (i=l,2) may change the conclusion. Especially if we chose Qi — Q2 ^ 180°, 
the two planets will in the aligned apsidal resonance , which is stable according to Fig.Qb. 
Similarly, HD37124 b and c could in apsidal resonance if suitable parameters Q are observed. 

Similar secular resonance may also happen due to nearly the same averaging precessing 
rate of the ascending nodes between the two planets in mutually inclined orbits. Since 
the inclination perturbations are isolated from the eccentricity ones in the linear secular 
perturbation theory, we will address this problem in a separate paper (Zhou and Sun, in 
preparation). 
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is supported by the Natural Science Foundation of China (No. 10233020), the Special Funds 
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Education of China for Ph.D candidate training(20020284011). 
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Fig. 1. — The approximation of cq in equation (3) 
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Fig. 2. — The critical case of HD12661 with the observed ai and q but with 02 = Oi?^ ~ 
1.76AU. (a) The phase plane 62 — Aw, the orbits in the diagram have the same constant 
D as that of eio = 0.35, 620 = 0.20. (b)Variation of 61,62, Aro with time for an orbit with 
Azuo = tt/2. 
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Fig. 3. — Libration region in 620 — Acuq plane defined by (19) and (20) with different Cio for a 
and q equal to the planetary system (a) HD 12661, and (b) 47 Uma. The black dots show the 
present configuration of the two planets in both systems, around which the boxes show the 
uncertainties of the orbital elements in Table 1 and 2. AtUa, A-oj^ are the lower integration 
limits in equations (23) and (26), respectively, and A-Ccj^, Aro^ are defined in equations (22) 
and (25), respectively. 
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Fig. 4. — Orbital diagrams of planet system HD12661 with initial values (a) eio = 0.1,Aci7o = 
0; (b)eio = 0.1,Aroo = tt 
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Fig. 5. — Contours of Aei, Ae2 defined by (14) in tfie initial 620 — Acuq plane with eio = 0.35. 
The dotted lines are the boundary of the libration region defined in (19) (20). 




Fig. 6. — Variations of the libration area ratios with a and q for the planet system HD12661. 
(a) the observed q ^ 1.46 is fixed, and (b) the observed a ^ 0.32 fixed. In both diagrams 
eio = 0.35 is fixed. The black dot in each plot shows the location of observed configuration. 
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Fig. 7. — Libration regions in the 620 — AciJo plane in the general three-body system for 
the HD12661 system. The solid curves show the boundary defined in (19)(20). The initial 
eccentricity is (a) eio = 0.10 and (b)eio = 0.35. 
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Fig. 8. — The LCE x(t) for (a) the 47Umas system with the observed orbital parameters 
listed in Table 2; (b) the HD12661 system with the observed orbital parameters listed in 
table 1; (c) same parameters as in (b) except 620 = 0.20, wiq = 0, ZU20 = O.Gtt; (d) same 
parameters as in (b) except 620 = 0.35, zuiq = 0, ZU20 — O.Ttt, the outer planet is escape at 
time t ^2.1 X lO^r. 




Fig. 9.— The LCE xe for HD12661 system. The initial eccentricity is(a) eio = 0.10, 
(b)eio = 0.35. The sohd curves show the boundary defined in (19) (20). The black dot in 
(b) shows the location of observed configuration with Aotq = 130.7°, around which the box 
show the uncertainties of the orbital elements. 




Fig. 10. — (a) The index of libration region and (b) the relative semi-major axis of two 
planets HD12661 systems under azimuthal acceleration f2 — —2 x 10~^AU^yr~^ . The 
initial eccentricity is eio = 0.35 and Mio,M2o are random chosen. The time evolution is 
t = 50000yr. 
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Fig. 11. — (a) The index of libration region and (b) the relative semi-major axis of two 
planets HD12661 systems under azimuthal acceleration /2 = —2 x lO^^AU^yr"^. The initial 
eccentricity is eio = 0.35 and Mio, M20 are random chosen. The evolution time is t = 
-SOOOOyr. 




Fig. 12. — Libration region in the initial 620 — ^-^0 plane defined by (19) and (20) for the two 
triple-planet systems (a)55 Cancri {zub — 99°,OTc = 61°, zu,^ = 201°); (b) Upsilon Andromedae 
{wb — 73^, Wc — 250^, ujd — 260°). Invisible boundary curves are out of the range of 620- The 
black dots show the present configuration of the two planets. Other orbital elements are 
taken from Fischer et al. (2003). 
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Fig. 13. — Libration region in the initial 620 — Awq plane defined by (19) and (20) for the 
planetary system (a) GJ876 {wb = 330°, = 333°),HD82943 {zub = 96°, zu^ = 138°); (b) 
HD37124 {zub = 60°, gj^ = 259°). The black dots show the present configuration of the two 
planets. Other orbital elements are taken from Fischer et al. (2003). 
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Table 1. Orbital Parameters^ for HD12661(1.07Mo) system 



Parameter 



HD12661 b 



HD12661 c 



planet mass Msini (M^up) 
period P(days) 
Tp (JD) (days) 
semi-major axis a (AU) 
eccentricity e 

argument of pericenter ui (deg) 



2.30 
263.6(1.2) 
2,449,941.9(6.2) 

0.82 
0.35(0.03) 
293.1(5.0) 



1.57 
1444.5(12.5) 
2,449,733.6(49.0) 
2.56 
0.20(0.04) 
162.4(18.5) 



*Data from Fischer et al. (2003) 



Table 2. Orbital Parameters^ for 47UMa(1.03Mo) system 



Parameter 



47UMa b 



47 UMa c 



planet mass Msini (Mj„p) 
Period P(days) 
Tp (JD) (days) 
semi-major axis a (AU) 
eccentricity e 

argument of pericenter ui (deg) 



2.54 
1089.0(2.9) 
2,450,356.0(33.6) 
2.09 
0.061(0.014) 
171.8(15.2) 



0.76 

2594(90) 
2,451,363.5(495.3) 
3.73 
0.005(0.115) 
127.0(55.8) 



''Data from Fischer et al. (2002) 
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Table 3. Extensions of Aotq in apsidal resonance for the observed systems 



Planet Pair 




^ ^1 
a = — 


ein 

620 


aligned Aotq 


anti-aligned Aroo 


Ups And b-c 


0.358 


0.0720 


0.037 


(-79.3°,79.3°) 


_c 


Ups And b-d 


0.181 


0.0228 


0.040 


(-47.0° ,47.0°) 




Ups And c-d 


0.507 


0.317 


1.08 


(-9.0°,9.0°) 




55 Cnc b-c 


4.15 


0.477 


0.073 




(96.8°,263.2°) 


55 Cnc b-d 


0.225 


0.021 


0.107 






55 Cnc c-d 


0.054 


0.044 


1.46 






GJ876 c-b 


0.296 


0.628 


2.70 




(144.0°,216.0°) 


47 UMa^ b-c 


3.34 


0.560 


12.2 


(-87.9°,87.9°) 




HD37124 b-c 


0.860 


0.184 


0.250 


(-69.8°,69.8°) 




HD12661 b-c 


1.46 


0.320 


1.75 


(-67.8° ,67.8°) 


(98.6°,261.4°) 


HD82943 c-b 


0.540 


0.628 


1.32 


(-59.7°,59.7<') 


(132.8°,227.2°) 


HD168443 b-c 


0.450 


0.103 


2.65 






HD38529 b-c 


0.061 


0.035 


0.806 






HD74156"^ b-c 


0.208 


0.080 


1.625 







■^Here - means no possible libration Aroo. 
'^Ba.ta from Fischer et al. (2002) 

''Data from California and Carnegie Planet Search (2003) 



